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Iron is the most important metal for modern industry and Sweden is by far the largest iron-producer in 
Europe, yet the genesis of Sweden's main iron-source, the 'Kiruna-type' apatite-iron-oxide ores, remains 
enigmatic. We show that magnetites from the largest central Swedish 'Kiruna-type' deposit at Grangesberg 
have 5 ls O values between —0.4 and +3.7%o, while the 1.90—1.88 Ga meta-volcanic host rocks have 5 18 0 
values between +4.9 and + 9%o. Over 90% of the magnetite data are consistent with direct precipitation from 
intermediate to felsic magmas or magmatic fluids at high-temperature (5 18 O mgt > +0.9%o, i.e. 
ortho-magmatic). A smaller group of magnetites (5 18 O mgt S +0.9%o), in turn, equilibrated with high-8 ls O, 
likely meteoric, hydrothermal fluids at low temperatures. The central Swedish 'Kiruna-type' ores thus 
formed dominantly through magmatic iron-oxide precipitation within a larger volcanic superstructure, 
while local hydrothermal activity resulted from low-temperature fluid circulation in the shallower parts of 
this system. 

Sweden is by far the biggest producer of iron in Europe 1,2 and iron ore has traditionally been sourced from 
two principal regions, the Kiruna-Malmberget province in northern Sweden and the Bergslagen region in 
central Sweden 3 (Fig. 1). The apatite-iron-oxide ores of the Grangesberg Mining District (GMD) represents 
the largest iron ore accumulation in the classic Bergslagen ore province 4 , with a past production of 156 Mt of ore, 
averaging 60% Fe and 0.81% P. Besides iron, phosphates such as fluorapatite and REE-silicates constitute an 
economically significant reserve of REE 5 . The apatite-iron-oxide deposits of the GMD and its northern con- 
tinuation are classified to be of 'Kiruna-type', whereas thousands of smaller, but still significant, iron-oxide 
deposits of banded-iron and skarn-iron types also occur in the Bergslagen province 4,6 . The GMD deposits 
themselves, are thus traditionally grouped together with those of the Kiruna-Malmberget region 6-8 (i.e., the type 
locality of the internationally renowned but genetically ambiguous 'Kiruna-type' ores). The origin of this ore-type 
has been a matter of substantial debate for over 100 years, and several fundamentally different modes of forma- 
tion have been suggested. These include direct magmatic segregation or crystallization, magmatic hydrothermal 
replacement, and hydrothermal precipitation in the sense of iron-oxide-copper-gold (IOCG-type) deposits 915 . 
Here we address the primary question, using oxygen isotopes, whether the GMD iron oxides formed through 
direct magmatic precipitation from either a magma or from high-temperature magmatic fluids (s 800°C), likely 
leading to intrusive-style, massive iron ores, or whether they formed through precipitation from aqueous hydro- 
thermal fluids at lower temperatures (e.g., £ 400°C), that would likely form disseminations and vein-type 
magnetites. 

Host rocks to the GMD ores comprise mainly intermediate to felsic meta-volcanic rocks that formed between 
~ 1 .90 and 1 .88 Ga in a subduction- or back- arc-related tectonic setting 4 . Both the apatite-iron-oxide ores and the 
meta-volcanic host rocks have been affected by later ductile deformation and the 1.85 to 1.80 Ga Svecokarelian 
low- to medium-grade metamorphism 4,5 . The exact timing of ore formation in the GMD is not known, but a 
major regionally-deformed granite body (> 1.85 Ga) occurs directly east of the main Grangesberg ore field and 
post-dates the ores (Fig. 1), because it carries characteristic metre-sized xenoliths of apatite-iron-oxide miner- 
alization. In most parts of the Bergslagen ore province, the host rocks to the ores are Si0 2 -rich dacites to rhyolites. 
In contrast, the host rocks to the GMD ores exhibit more intermediate compositions with a dominance of 
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Figure 1 | (A). Overview map of Fennoscandia with the Grangesberg mining district (GMD), the Bergslagen province, and the Kiruna-Malmberget 
mining district indicated. (B). Geological map over the main ore zone in the GMD. (C). Vertical section (line X-Y in panel B) through the main ore body 
of the GMD. The ore zone extends downwards at a 70-80" dip to the SE. Black horizontal lines are adits. Modified from Geological Survey of Sweden 
(SGU) datasets and A. Hallberg, SGU. 
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Figure 2 | Chondrite-normalised REE diagram of host rocks and iron-oxide ores from drill cores 717 and 690. Chondrite values for normalisation after 
Sun & McDonough 33 . The iron-oxide ores and the moderately altered host rocks share similar REE concentrations and patterns with the regional bedrock 
and recent subduction zone volcanic rocks. The hydrous altered host rocks (marked with hatched pattern), in turn, are enriched in REE relative to 
moderately altered host rocks and recent subduction zone volcanic rocks 36,37 . 
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Table 1 Grdngesberg oxygen isotope data, including sample position along the drillcore 
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andesites and dacites. The major ore body itself consists of several 
lenses that occur along a planar NNE-trending zone with a moderate 
to steep eastward dip. The principal orientation of the major ore body 
is stratiform, i.e. it is parallel to the internal contacts within the meta- 
volcanic host rocks that themselves are sub-parallel to the main 



tectonic trend in the region. Magnetite occurs either in massive ore 
lenses or as veins, bands or disseminations (hereafter referred to as 
VeDi-magnetites). Within the central part of the GMD ore field, a 
number of dykes and sills exhibit cross-cutting relationships to rock 
units of similar composition, on which basis we interpret the 
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Figure 3 | Parts (A) and (B) show the oxygen isotope data for two drill cores (Numbers 690 and 717) that traverse the main ore zone at Grangesberg 
between 570 and 670 m below the surface (see Fig. 1 ) . Shown are the oxygen isotope compositions of the host rocks, quartz separates, massive magnetites 
and VeDi-magnetites, the latter including magnetite from hematite ore. All oxygen data are reported in standard 5 ls O-notation relative to SMOW after 
Hoefs 34 . The S 18 0 ranges for the mantle and arc-andesites are after Bindemann 18 and Taylor". Range of igneous magnetites after Taylor 35 . 



emplacement environment as "shallow sub-volcanic" in character. 
Hydrous alteration is evident in the direct host rocks to the ore in 
the form of disseminated and sometimes discrete biotite, chlorite 
or variably amphibole- and fluorapatite-rich assemblages. Massive 
magnetite ores seem unaffected by this alteration or by the regional 
metamorphism. However, the hydrous alteration assemblages are 
always spatially related with the massive ore and are absent beyond 
the main ore zone, implying a genetic relationship between the two. 

Results 

Samples were collected from surface exposures and from three shal- 
low plunging (s 20°) drill cores that transect the steeply dipping ore 



body between —570 and —670 m (Fig. 1). The GMD apatite-iron- 
oxide ore is either magnetite- or hematite-dominated, and is variably 
enriched in fluorapatite and accessory silicate phases (amphiboles, 
biotite, chlorite sensu lato) and REE-minerals. Most ore-types show 
some degree of banding, defined by aggregates of mostly fine-grained 
fluorapatite and minor amounts of silicates. The main GMD ore 
field, locally known as the 'Export Field', is estimated to consist of 
about 80% magnetite and 20% hematite ore, with a concentration of 
the latter near the structural footwall. 

In addition to fluorapatite, REEs are hosted in associated mon- 
azite-(Ce), allanite-(Ce), xenotime-(Y), LREE-bearing epidote, and 
REE fluorocarbonates. Although the GMD apatite-iron-oxide ores 
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Table 2 | Calculated oxygen isotopic compositions for magmas and fluids in equilibrium with magnetite samples 

Theoretical magma Theoretical magma Theoretical aqueous ore 

5 ls O (%») for 5 ,8 0 (%o) for fluid 5 ,s O (%■>) for 
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y = in equilibrium with magma/magmatic water ; X= not in equilibrium with common magmatic values. 



are typically enriched in Th, U, La, Ce, Nd, P, Fe, Sm, Tb, Y, Tm and 
Yb, they are depleted in K, Ba, Sr, Zr and Ti relative to average 
continental crust. The ores themselves, and particularly the apat- 
ite-rich samples and the hydrous host rock alteration assemblages, 
display negative Eu-anomalies and flat HREE patterns. The hydrous 
altered host rocks and the apatite-rich assemblages associated with 
the ore show the same REE enrichment patterns as moderately 
altered host rocks, but are elevated by several orders of magnitude 
when compared with the meta- volcanic rocks in the region and with 
typical subduction zone andesites and rhyolites (Fig. 2). 

To distinguish between magmatic and magmatic-hydrothermal 
origins of magnetite, which together we term "ortho-magmatic", is 
not possible with our data. However, a low-temperature hydro- 
thermal origin 400°C) versus a high-T ortho-magmatic origin 
can be discerned. Oxygen isotope ratios have been measured on 63 
samples of meta-volcanic host rocks, hydrous alteration-assem- 
blages, quartz mineral separates, as well as magnetite from massive 
ore, veins, and disseminations, and from hematite-dominated ore. 
All data are reported in delta notation relative to SMOW (Table 1). 
The host rock 8 18 0 values, including whole rocks and quartz mineral 
separates from both outcrops and drill cores, range from +4.9 to 
+ 10.3%o, i.e. within the spectrum of most normal igneous rock 
compositions. The massive iron-ores from the main ore body yield 
5 ls O values, between +0.9 and +3.0%o [n = 16] (Fig. 3). The S ls O 
values of all GMD magnetite samples, including vein and dissemi- 
nated magnetite types [n = 7] as well as magnetite inclusions in 
hematite ore [n = 2] range from —0.4 to +3.7%o. The origin of 
the different magnetite types can be assessed using O-isotope frac- 
tionation between magnetites and an andesitic parent magma at 
magmatic temperatures (> 900°C) (A magnetite _ andesite = -4.0%o) 
or with magmatic water (A magnetite _ water = — 5.2%o at 800°C) 16,17 . 
Magnetites with 8 18 0 £ +1.7%o yield equilibrium magma values 



between +5.7 and +7.7%o (Table 2), i.e. within the range of subduc- 
tion zone andesite magmas ( + 5.7 to +8.5%o) 1819 . Using the frac- 
tionation factor for magnetite from a dacite magma 
(Am. gnetite— dacite — 4.3%o) 17 , the theoretical magma values would 
shift upwards by +0.3%o, moving them further into the intermediate 
magmatic range that is characteristic for continental subduction 
zone rocks 1819 (Figs. 3, 4, Table 2). These results are compatible with 
our andesite to dacite host rock data that range between +4.9 and 
+ 9.0%o (Table 1). In turn, magnetite samples with 5 ls O values s 
+ 0.9%o satisfy equilibrium conditions with high-temperature mag- 
matic waters of +6.1 to +8.9%o (Table 2), which corresponds to the 
range of common magmatic waters (6-8%o) 20 . Thus, all GMD mag- 
netites with 8 18 0 values s +0.9%o (> 90% of our samples) are 
consistent with an 'ortho-magmatic' mode of formation because they 
satisfy high-temperature magmatic equilibrium conditions (pale 
pink box in Fig. 4). Notably, this field hosts all GMD massive mag- 
netite ore samples, as well as several VeDi samples. Magnetites with 
8 18 0 values below +0.9%o, in turn, cannot be in equilibrium with a 
regular magma or with a high-temperature magmatic fluid, but must 
have formed from a high-8 18 0 hydrothermal fluid at more moderate 
temperatures. Low-temperature processes increase mineral-fluid 
fractionation, i.e. fractionation factors become more negative (e.g. 
Amagnctite-water — ~ 7.6%o at T < 400°C) 16 , implying oxygen isotope 
exchange at low temperatures for the magnetites £ 0.9%o. A high- 
8 18 0 hydrothermal fluid from the uppermost range of meteoric 
fluids 21 , or a mix of remnant magmatic water and a meteoric influx, 
can explain the remaining VeDi magnetites (n = 2), while satisfying 
the low-temperature fractionation conditions (Fig. 4). 

Employing a similar approach on host rock quartz, a mineral 
that is assumed to be resistant to oxygen isotope resetting once 
formed, should allow us to test the results derived from the mag- 
netite data. Quartzes have 8 ls O values between +5.8 and +10.3%o 
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Figure 4 | Magnetite 5 is O values from GMD compared to other volcanically-hosted iron ore deposits. For reference, magnetites in equilibrium with 
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evolved rhyolite with a S ls O of 10%o (the demarcation between I-type (<10%o) and S-type (>10%o) magmas) are shown. The GMD magnetites plot 
dominantly above the +0.9%o demarcation and in the field of 'ortho-magmatic' magnetites after Taylor 35 ), and satisfy equilibrium with magma or 
magmatic fluids at magmatic temperatures (~800-1000°C). A small fraction of the GMD data (n = 2), however, is more consistent with formation from 
a low-temperature fluid regime. The cut-off point for this is calculated to be +0.9%o in magnetites, because fractionation factors determine that samples 
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(n = 11, Table 3). Taking A qtz _ magma = +0.9%o and +0.8%o 17 as 
fractionation between quartz and andesite and quartz and dacite 
magmas respectively, 5 ls O melt values of +4.9 to +9.5%o are 
derived (n = 22, Table 3). Seven out of 11 quartz 8 ls O values 
correspond to calculated magmatic values for andesites and 
dacites with S 18 0 s + 5.7%o 18 . A calculated fractionation temper- 
ature for a co-existing quartz and magnetite pair for the GMD 
deposit that appears in textural and oxygen-isotope equilibrium 
yields a temperature of 907 ± 53°C, consistent with temperatures 
expected for oxide formation from intermediate to evolved (hyd- 
rous) arc magmas and consistent with phase relations and experi- 
mental constraints 9 ' 10 - 13 ' 22 . Quartz with 5 ls O < +6.2%o (n = 4) 
would be in equilibrium with a magmatic fluid at temperatures 
between 800°C and 900°C (A quartz _ water = +0.3 to -0.1%o), 
which implies that these quartzes are of ortho-magmatic origin 
too (following our initial definition). 

Discussion 

The relatively wide range of 5 ls O in the GMD host rocks ( + 4.9 to 
+ 9.0%o, Table 1) suggests that some host rock and VeDi-samples 
reflect low- and high-temperature alteration processes (e.g. silicate 
whole rocks with 8 18 0 > 8.5%o and < 5.7%o respectively). Massive 
magnetite ores, in turn, show textural characteristics and 8 ls O values 
that support precipitation from a magma. Those massive magnetites 
are chemically and mechanically refractory and appear little affected 



by later overprint (i.e., they were resistant to subsequent exchange 
with fluids or there were insufficient fluids to exchange with). In 
addition, the temperature and oxygen isotope data derived for the 
refractory GMD quartz and magnetite samples fall in the magmatic 
range also and overlap with those presented in Nystrom et al. 20 for 
the iconic El Laco apatite-iron-oxide deposit in Chile (—700° - 
800°C). Discussion is, however, ongoing regarding the genesis of 
the El Laco deposit, and the possibility of interaction with a high- 
8 18 0 hydrothermal fluid exists because the data for almost all textural 
magnetite types at El Laco average around +4%o 23 . These authors 23 
argue that a 'natural' spread between "magmatic" and "low-temper- 
ature hydrothermal" signals would be expected for a truly "mag- 
matic" formation because a large volcano would logically develop 
an associated hydrothermal system at shallower levels. Although 
maybe not observed at El Laco, the GMD data show a small group 
of low-temperature magnetites, making our proposed dominantly 
ortho-magmatic formation scenario rather probable. The potential 
wider significance of this result is further underlined by the fact that 
our GMD data encompass the available 8 ls O-data from Kiruna 20 , the 
largest apatite-iron-oxide deposit in Europe, suggesting very similar 
modes of formation for the GMD and Kiruna deposits (Fig. 4). 

On the basis of the data presented, the massive apatite-iron-oxide 
ores at GMD formed in equilibrium with mafic to felsic magmas at 
magmatic temperatures (i.e., probably with subduction-type ande- 
sites to dacites, which form the host rocks to the ores in the region). 
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The hydrous alteration assemblages, in turn, have highly variable 
8 18 0 values, but their REE patterns are very similar to that of the 
ores. In combination with their exclusive occurrence close to the ore 
zone, this observation suggests a genetic relationship between iron- 
oxide ores and the REE -enrichment in hydrous alteration assem- 
blages 2426 . Medium- to low-temperature hydrothermal alteration 
associated with ore formation and emplacement is consistent with 
a scenario of high-temperature oxide-rich magmas intruding an 
upper crustal volcano-sedimentary pile where they initiate fluid cir- 
culation at more moderate temperatures. At GMD, volcaniclastic 
host rocks are locally extensively brecciated as well as crosscut by 
dykes and other sheet intrusions, which in combination with our 
REE and oxygen isotope data supports a shallow-level sub-volcanic 
environment. The VeDi-type magnetites and hydrous alteration 
assemblages formed in proximity to the massive ore bodies, and a 
degree of overlap between these regimes (i.e., high-temperature 
ortho-magmatic versus low-temperature hydrothermal), is expected 
for a sub-volcanic mode of formation 23 . Our preferred interpretation 
is therefore that of a now deformed, originally large and possibly of 
caldera-type 4 , volcanic complex situated in a subduction zone 
regime. Ascent, intrusion, and crystallization of oxide-rich melts to 
shallow levels within the volcanic plumbing system produced dom- 
inantly ortho-magmatic magnetite, while low-temperature hydro- 
thermal activity near the surface occured to some extent. This late 
hydrothermal activity likely involved exsolved magmatic volatiles as 
well as external meteoric waters that probably became available to the 
host system at shallower levels and at progressively lower tempera- 
tures 23,27 . This realisation helps us to better understand the formation 
of Europe's most important iron source, the 'Kiruna-type' apatite- 
iron oxide ores of central and northern Sweden. 

Methods 

REE analysis. A total of 44 drill core samples including both, moderately to intensely 
and hydrous altered host rocks (9 and 15 samples, respectively) and sections from the 
iron oxide ore (20 samples) were analysed for Rare Earth Elements (REE) at Acme 
Analytical Labs Ltd in Vancouver, Canada. Trace and rare earth elements were 
analysed by inductively coupled plasma-mass spectrometry (ICP-MS) after 
preparation by multi-acid digestion. Loss on ignition (i.e. volatile content) was 
determined by igniting a sample split and then measuring weight loss. Data quality 
was monitored using a set of certified, internal reference materials. Detailed 
information on uncertainties and reproducibility can be found at http:// 
acmelab.com/services/method-descriptions/soil-till-and-sediment/. 

Oxygen Isotopes. All O-isotope data presented in this paper were produced at the 
University of Cape Town (UCT). For oxygen isotopes, both conventional and laser 
fiuorination methods were used. Some quartz, and all magnetite, separates were 
analysed using a conventional silicate line (described in Harris & Ashwal 28 ). 
Approximately 20 mg (magnetite) and 10 mg (quartz) of sample was reacted with 
CIF3, and the liberated 0 2 converted to C0 2 using a hot platinized carbon rod. 
Magnetite was reacted overnight at 600°C and all other minerals were reacted for 
4 hours at 550°C. Some quartz and magnetite separates were analysed using the laser 
fiuorination analytical process described in Harris & Vogeli 29 . Each sample was 
reacted in the presence of approximately 10 kPa BrF 5 and the purified 0 2 was 
collected onto a 5A molecular sieve contained in a glass storage bottle. 

All isotope ratios were measured off-line using a Finnigan Delta XP mass 
spectrometer in dual-inlet mode. All data are reported in 5 notation where 8 18 0 — 
(Rsampie/Rstandard ~ I)*1000, and R = the measured ratio (i.e. ls O/ 16 0 or D/H). 
Duplicate splits of the quartz standard (NBS28) run with each batch of eight samples 
were used to convert the raw data to the SMOW scale using the 8 ls O value of 9.6%o for 
NBS28 recommended by Coplen et al. 30 . During the course of this work, 8 analyses of 
NBS28 gave a 2a error of 0.2. The O-isotope ratios of samples analysed using laser 
fiuorination were measured on 0 2 gas. Measured values of our internal standard 
MON GT 31 were used to normalise raw data and correct for drift in the reference gas. 
The average difference in 8 is O values of duplicates of MON GT analysed during this 
study was 0.1%o, and corresponds to a 2a value of 0.2%o. MON GT was recalibrated 
against the UWG-2 garnet standard of Valley et al. 32 using the current laser system, 
and has a revised 8 18 0 value of 5.4%o, assuming a 8 18 0 value of 5.8%o for UWG2. 
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